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Calculations of the pairing interaction in multi-band models of the Fe superconductors show that 
it is attractive in both the Aig (s-wave) and -Big (d-wave) channels. This raises the possibility that 
these materials may have collective excitonic modes. Here, assuming an s-wave groundstate, we 
investigate the d-wave collective excitonic mode and its coupling to the Raman scattering. 



PACS numbers: 

Both RPA fluctuation exchange^'^ and numerical func- 
tional renormalization group calculations'^ find that s- 
wave (Aig) and d-wave (Big) instabilities can occur in 
multi-band models of the Fe-pnictides. Typically the 
pairing strength or Tc of the s-wave state is found to be 
greater than that of the d-wave state. Nevertheless, there 
are parameter ranges where the two states lay relatively 
close, raising the possibility that there could be a d-wave 
{Big) collective excitonic mode.^'^ In a fully gapped su- 
perconductor, this mode, consisting of two quasi-particles 
in a "Cooper pair" d-wave state, lays below 2A for zero 
center of mass momentum and would appear as a sharp 
peak. For the Fe superconductors, calculations suggest 
that the s-wave gap in the groundstate is anisotropic 
and may even have nodes. In this case, the collective 
state could be damped and appear as a broad resonance 
in the two quasi-particle d-wave channel. Experimen- 
tally such a mode with L = 2 could be excited from an 
s-wave superconducting state by Raman scattering. 
The observation of such a resonance would be the first for 
such a collective state. Early microwave measurements 
found indications of a precursor response just below 2A 
in Pb and it was initially thought that this might be a 
p-wave collective mode. However, further studies deter- 
mined that this was an artifact.^ Here we explore the 
possibility of the existence of a d-wave collective mode in 
the Fc-pnictidcs and examine how it could be detected 
by Raman scattering. 

We will begin by considering the simple case illustrated 
in Fig. la. Here we imagine that the two hole Fermi sur- 
faces around the F point of the 1 Fe/cell Brillouin zone 
have been collapsed into one a Fermi sheet and Pi and 
/?2 represent the two electron Fermi sheets. Suppose the 
s-wave (Aig) part of the pairing interaction connecting 
the a and /3 Fermi sheets is denoted by V with a cut- 
off LjQ on I Ek I around the Fermi surfaces. There will 
of course also be (3i — [32 contributions to the s— wave 
pair scattering as well as intra a, [3i and (32 terms which 
we neglect in this simple model. Then the Aig pairing 
instability represented by the diagrams in Fig. lb is de- 
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FIG. 1: (a) The Fermi surfaces in the 1 Fe/cell Brillouin zone. 
Here we have combined the ai and Q2 hole Fermi surfaces into 
one a sheet. For the simple example that we will consider 
the pairing interaction for the A\g (s— wave) ground state 
is desribed by a strength V for scattering pairs between the 
a and /3i and {32 Fermi surfaces. The d— wave part of the 
interaction U scatters d— wave pairs between /3i and (32. (b) 
Diagrams for the A\g pairing instability. 
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with Gf3^a{k,iiOn) — {ii^n — e0{k))~^ the single particle 
Green's function on the (3 or a Fermi surfaces, V the a-/3 
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pairing interaction, and (f)fj,a the gap function amplitude 
on the (3 and a Fermi surfaces, respectively. From Eq. 
(1), the transition temperature is determined by 
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w ith the s- wave coupling strength Ag = 
sj2Ni3{Q)Na{Q)V. The ratio of the gap amplitudes at 
is 
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In the superconducting state, the dominant d-wave 
scattering between the quasi-particles occurs between the 
/3i and (52 Fermi surfaces. We will parameterize the d- 
wave part of this interaction by the separable form 

r<^(fc, k') = -gl {6)Ugl {B') - {6)Ugl {9') (4) 

for |£fc| and |efc'| less than a cut-off frequency ujq. Here, 
gfjXO) depends upon the angle of k on the /3i-Fcrmi sur- 
face measured from the kx axis and 5^ {9) = —g'^ {9 + 

Similarly to our treatment of the s— wave pairing chan- 
nel, here for simplicity wo will neglect /? — a contributions 
to the d— wave channel. The additional contributions to 
both the s— wave and d— wave pairing interactions basi- 
cally only change the strengths of the effective s— and 
d-wave pairing interactions and A^ which we take as 
parameters in the following. 

If the system were to remain in the normal state, su- 
percooled below the ,s-wavc pairing instability, it would 
become unstable to pairing in the d-wave channel when 
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Normalizing the angular average of gj{9f to unity 
around the ^\ Fermi surface, the d-wave transition tem- 
perature is 
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with the d-wave coupHng strength A^ = N^{0)U. Fluc- 
tuation exchange calculations^'^ and numerical renormal- 
ization group studies for models of the Fe superconduc- 
tors find that the coupling strength Ad can be comparable 
to the coupling strength in the s-wave channel As, raising 



the possibility that one may find a d-wave "Cooper Pair" 
collective mode in the s-wave superconducting state. 

The homogeneous Bethc-Salpeter equation for a col- 
lective d-wave mode is illustrated in Fig. 2a. Here 
k = (k, iw„) and q = (q, zwm)- The single- and double- 
arrow lines denote the single particle Green's function 
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and Gor'kov's anomalous Green's function 
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in the s-wave superconducting state, respectively. The 
Bethe-Salpeter equation for the d-wave collective mode 
is 
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X {G{k' + q)G{k') + F{k' + q)F{k')) (8) 

For the separable interaction given by Eq. (4) one obtains 
an equation which determines the energy of the collective 
d-wave mode 



k'n' 



X [G{k' + q)G{k') + F{k' + q)F{k')] . (9) 

After the Matsubara sum is evaluated and T ^ 0, we set 
q = and analytically continue icUm + iS to give 
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with E = ^e'^ + A'^{9) and w is assumed to have a 

small positive imaginary part. Here A^(^) is the s-wave 
groundstate gap on the /3 Fermi surface which we will 
set equal to Aog^{9). The integral over e is done in the 
usual way 
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and extending the range of integration for the first term 
to plus and minus infinity gives 
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FIG. 2: (a) Diagrams for the d— wave collective mode. Here 
the thick single arrow line represents the single particle Greens 
function G(p,ia;„) and the double arrow line the Gor'kov 
anomalous Greens function F{p,iu!n) in the Aig state, (b) 
The lowest order U contribution to the Raman scattering from 
interband pair interactions, scattering pairs between /3i and 

02. 



The collective d-wave mode at g = has a frequency 
and damping given by 
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where the bracket implies an angular average. Here the 
tilde s-wave coupling strength is 
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Depending on the difference in coupling strengths and the 
anisotropy of the s-wave gap on the (3 Fermi surfaces, one 
will have a sharp mode or a resonance. 

A first order contribution of the interaction vertex 
to the Raman scattering is illustrated in Fig. 2b. There 
are four arrangements of the G and F propagators and 
one can go from fii to /32 or (32 to pi . Finally adding the 
spin sum, the first order contribution of F'' to the Raman 
susceptibility is 



Ax(*^™) = {/(47^,Gi^.9^J(%,Gi^7&) (15) 



with 
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FIG. 3: Plots of the Raman response for different gap 
anisotropics r in the absence of collective mode effects (black) 
and with collective contribution included (red) for As = 1 and 
Xd = 0.5. Here the parameter a has been set to 1 and a small 
damping term has been added. Note the changes in scale of 
the r = plot. 



Here q = {q, iujm) and we are interested 'm q — Q. Evalu- 
ating the Matsubara sum, we have 

(47^GFg^) - 7V^(0)Ao(7^(%^(^)<?Me)P(^,e)) (17) 



and 
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Here again the bracket implies an angular average. 

Symmetry considerations can be used to determine the 
collective mode contributions from the interplay of domi- 
nant and sub-dominant pair interactions and polarization 
geometries. From Eqs. 11, 12, 17 and 18, one can see 
from symmetry that the collective mode contribution to 
the Raman vertex will vanish unless g'^^jsg"^ transforms 

as Aig for tetragonal D^^ symmetry. As used here, since 
g^g^ transforms as one of the d— wave representation, a 
collective mode will appear only for crossed polarization 
incoming and scattered polarization geometries. Specifi- 
cally, if one considers V to be in the (1^2 _y2 channel, the 
collective mode will appear for Big orientations only. 
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FIG. 4: Same as Fig. 3 but for = 1 and Ad = 0.8. 



As previously discussed, multiple scattering between 
the f3i and P2 Fermi surfaces leads to a collective d-wave 
Cooper pair state. Its contribution to the Raman scat- 
tering is obtained by replacing U in Eq. (15) with 
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Proceeding as before, we find that 
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The lowest order contribution to the /3-Fermi surface 
Raman scattering is given by^'' 
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If we set {'^pg'ljY = '^{ifj) ™ Eq. (21) and (22) we have 
the following expression for the Raman scattering 
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Here we have also set {gi{0)'^) = 1. 




Figure 3 shows plots of Imx/N;3(0)(7|) for Apie) = 
Ao(l + r cos26')/(l + r) for various values of r and cou- 
pling strengths (As = IjA^ = 0.5). The results indicate 
that the collective mode removes spectral weight from the 
higher energy portion of the response and adds weight at 
the mode position, determined by l/A^ — l/Ag. For the 
isotropic gap (r — 0), an essentially undamped mode ap- 
pears below 2A, capturing most of the spectral weight 
from the bare contribution (the first term in Eq. 23). 
For a gap anisotropy r = 0.6 that still preserves a true 
gap, a well-defined collective mode appears at frequen- 
cies slightly below 2Amm, removing the log singularity of 
the bare response at 2Amax and shifting spectral weight 
to lower energies. For a gap with nodes (r = 1.4), the 
collective mode is damped by the finite particle-hole con- 
tinuum at all Raman energies, and re-allocates spectral 
weight to low energies near | 2Ami„ |. 

For a stronger d— wave couplings Xd = 0.8, appropriate 
for the case of near-degenerate pair interactions, the col- 
lective mode is pulled further away from the continuum 
contribution, as shown in Fig. 4. For the isotropic gap 
case, the mode drops to lower energies and has a smaller 
residue. This is the case for r = 0.6 as well, where a 
well-defined collective mode pulls out of the continuum 
below 2Amin- For the nodeful case r = 1.4, the col- 
lective mode again drops to lower energies but remains 
damped. It however changes the low-frequency behavior 
of the Raman response considerably. 

For the Fe-pnictides, contributions to the Raman re- 
sponse comes from each of the a and /3i_2 bands, plus 
mixing terms. As the d-wave collective mode contribu- 
tion arises only from multiple scattering among the /3 
bands, the contribution from the a bands will be deter- 
mined from simple non-interacting considerations given 
in Ref. 10. For the case of degenerate pi, 2 bands hav- 
ing energy gaps A,3^^/3^{9) = Ao(l ± rcos26')/(l + r), 
respectively, the mixing terms vanish by symmetry for 
Big polarizations, and the collective mode contribution 
is twice Eq. 23, shown in Figs. 3 and 4. 

Since the collective mode contribution is predicted to 
appear in crossed initial and scattered photon polariza- 
tion orientations, this Big particle-particle mode will not 
be coupled to long-range Coulomb forces, which nomi- 
nally push Aig (parallel polarizations) collective mode 
contributions up in energy to the plasma particle-hole 
state. This may allow for this mode contribution to ap- 
pear at low energies distinctly separate from the contin- 
uum in Raman experiments on the pnictides, unlike the 
situation in conventional superconductors.^ 

Therefore the detailed lineshape of the electronic Ra- 
man continuum may be determined from the interplay of 
anisotropics of the Raman vertices and the structure of 
the pairing interaction using symmetry arguments, as ap- 
plied successfully in the cuprates^''''^^. What may distin- 
guish the Fe pnictides from the cuprates may be the pres- 
ence of pairing channels having almost equal strength, 
as indicated in recent spin-fluctuation and RG consider- 
ations. This paper has shown that such a circumstance 
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will result in collective mode contribution which will have 
an unique polarization signature in the Raman spectrum. 
This could open a window into the determination of the 
pairing structure of the pnictides and provide important 
clues to the pairing mechanism. 
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